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INELASTIC SCATTERING CALCULATIONS WITH PROJECTED HARTREE-FOCK 
WAVE FUNCTIONS - COUPLED CHANNEL TREATMENT 
by Robert J. Ascuitto, * Richard C. Braley, and William F. Ford 

Lewis Research Center 

SUMMARY 

20 24 

The inelastic scattering of protons by Ne and Mg is investigated using a tech- 
nique (coupled channels) that explicitly includes the effect of the low-lying nuclear states. 
These states are represented by many-nucleon wave functions obtained by extracting 
the appropriate rotational component from an intrinsic Hartree-Fock wave function. 

The cross sections obtained in this way are compared with others obtained in the usual 
sort of calculation (distorted wave Born approximation), in which the effects of all 
nuclear states are lumped into the optical potential. The results show that, although 
angular distributions for the first level can be calculated correctly using the simpler 
method, those for the higher levels are given accurately only in the coupled channel 
approximation, and not always then. The variation between the two methods is different 
for each nucleus, and the reasons for this are discussed and shown to support the valid- 
ity of the Hartree-Fock model. Finally, the lack of agreement between the experimen- 
tal and theoretical cross-section magnitudes is discussed, and it is suggested that the 
problem is not with the nuclear wave functions, but apparently (and rather unexpectedly) 
with either the two-body force or the reaction mechanism. 


INTRODUCTION 

In recent years, a great deal of interest has developed in understanding inelastic 
nucleon scattering in terms of the correlations responsible for the nuclear transition. 
Until very recently (ref. 1), calculations have been restricted to nuclei, which are well 
described by the "shell model" and "vibrational model" (refs. 2 and 3). In both de- 
scriptions, relatively simple (inert-core) nucleonic configurations are considered 
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dominant components in the wave functions, thereby simplifying the calculation of form 
factors needed to describe the inelastic scattering. 

A wealth of experimental data also exists for scattering from strongly collective 
nuclei, since these have the largest inelastic cross sections. Unfortunately, the en- 
hanced transition strength in these nuclei is an indication of strongly correlated motion 
among many nucleons, so that a description in terms of only a few active (extracore) 
nucleons is likely to fail. We therefore use a microscopic picture. By this we mean 
that the nucleus is described in terms of the interactions between all the nucleons. 
Furthermore, the scattering will be described microscopically in the sense that the 
interaction between projectile and target shall be a sum of two-body interactions between 
projectile and target nucleons. 

In an effort to examine how the participation of the "core nucleons" affects transi- 
tion strengths, Love and Satchler introduced a core polarization model (ref. 4). This 
procedure essentially enhances the usual transition amplitude by addition of a term that 
simulates the activity of the core by means of the collective model. A more complete 
microscopic description of the core participation in inelastic scattering has been given 
by Braley and Ford (ref. 1). Core polarization effects were included explicitly by 
allowing all of the deformed orbitals to vary in a Hartree-Fock (HF) calculation. The 
wave functions obtained from that calculation were subsequently used in a distorted 

wave Born approximation (DWBA) study of inelastic proton scattering to the first 2 + 

20 

state in neon-20 ( Ne). It was found that the major effect of core polarization for the 
0+ - 2| transition was to increase the magnitude of the cross section. 

In contrast to the microscopic model, the macroscopic collective model absorbs all 
knowledge of the underlying correlations by parameterizing this information in terms of 
surface deformations which can be related to various transitions of the nucleus (ref. 5). 

Although this approach has been extremely successful in reproducing experimental 
results, it leaves one with a very limited knowledge of the nucleus and the correspond- 
ing inelastic scattering process. In the (microscopic) HF approach, collectivity or de- 
formed shapes are explained as a consequence of the self-consistent field generated by 
the interacting nucleons Naturally, this is done in some approximation that neglects 
certain correlations. However, it has the appealing feature that deformation parameters 
may be retrieved while, at the same time, the nuclear state is described in terms of 
the individual nucleonic motion. Of course, the validity of the HF method rests largely 
on comparison with experimental data. Nevertheless, one might also expect to gain 
some insight into the inadequacies of the microscopic model by examining the interplay 
between it and the macroscopic model. 

Recent analysis of proton inelastic scattering on certain s-d shell nuclei, using the 
macroscopic model, showed evidence of strong hexadecapole shape deformation in 
addition to the usual quadrupole shape deformation (ref. 6). To reproduce these 
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moments, or the corresponding multipole transition strengths needed in the inelastic 
scattering, would be a strict test of the Hartree-Fock scheme. In the present calcula- 
tion, we have tried to see how well we could reproduce these experimental results using 
the Hartree-Fock method to generate the intrinsic states 0^. The ground-state rota- 
tional band is obtained by projecting states of good angular momentum from the intrinsic 
state. The s-d shell nuclei have the advantage that the number of nucleons is suffi- 
ciently small so that HF calculations can be carried out with all particles active (i.e. , 
no core). This allows us to study problems such as basis truncation in a convenient way 
since all particles are treated on an equal footing. 

A coupled channel (CC) calculation of proton inelastic scattering is carried out for 
neon-20 (^Ne) and magnesium-24 (^Mg). In each case, at least the 0 + , 2 + , and 4 + 
members of the ground-state band were included explicitly in the coupled equations. 

The CC results are then compared with results of the DWBA. In an attempt to make 
the comparison with DWBA meaningful a special effort is made to choose DWBA optical 
model parameters that yield the same elastic scattering as the CC study. This proce- 
dure should remove doubts that would normally be raised because of differences in the 
elastic scattering that would manifest themselves in the inelastic cross sections by 
means of the distorted waves. 

NUCLEAR STRUCTURE INFORMATION AND CALCULATIONAL PROCEDURE 

Structure 

20 24 

The wave functions for the states of interest in Ne and Mg are obtained by pro- 
jecting states of good angular momentum from an axially symmetric intrinsic HF state 
0 K ; that is, 

^JM P MK ( la ) 


where (ref. 7) 


I^MK _ r f&Sl (S7)*R (O) (lb) 

8 ir 2 J 

In the standard approach one generates a rotational band by projection from a single 
Slater determinant, 0^ 
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Most of the predicted nuclear properties that will be studied involve reduced matrix 
elements of sums of single particle operators, * that is, 


n = 


A 

^(x n ) 

n=l 



(4) 


When this is the case, the reduced matrix elements may be expressed as (ref. 8) 

( f 1 i 0 1 1 i) jm =^Sj(if|ab) (blifilia) JM (5) 

ab 

where a and b label single-particle states and 


Sj(if lab) =£ua a J'JJ f ; Jji b ) ( (6) 

J' 

Jl J' 

The reduced overlap integral ( involves initial and final nuclear states <p^ and 

(p^ with the single -particle states <p & and 0 b projected out. When the deformed HF 
method with projection (ref. 7) is used to obtain the initial and final nuclear states, the 
reduced overlap integral may be expressed as (ref. 1) 



XX' MM' 


x <J c M\j b m' J f K’> 


2J, 


-/ 


^M’M^ Me) @X'X ^ ) ® ^ 


where 

Me) - <^ K 'i e " ieJy !^K> 

P x . x (0) = (B ) AX . 

B xx ,(0) =<x|e _i0J y|x'> 

■^Expressions for the operators used in the present calculations may be found in 
ref. 5. 
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The calculation of the projected energy spectrum is more complicated since the 
nuclear Hamiltonian 3C contains two-body operators in addition to the one -body kinetic 
energy terms. One must evaluate 

E j = (^kI 3CP kkI^k^ // ^kI p kk1^k^ 

There are several methods that may be used to carry out the energy projection (ref. 7). 
The most efficient of these is one in which the rotations are carried out explicitly. 

This technique is discussed in detail elsewhere (ref. 9). 


Scattering 


Inelastic scattering based on microscopic nuclear descriptions lias been considered 
by several authors. The coupled- channel approach to inelastic scattering has been dis- 
cussed by Glendenning (ref. 2). This treatment involves the solution of a set of differ- 
ential equations subject to the boundary condition that there be outgoing waves in all 
channels c but incoming waves only in the elastic channel Cq: 



Hi + D1 


u S P t« - E < 


u 


n I 

cc 


o W = X < 

u c Vc 




c c, 


(r) 


( 8 ) 


where E is the channel energy and the functions 0 represent a vector coupling of the 
spin-orbit functions of relative motion (Y^ g j) and the nuclear state functions (0 j^yj) to a 
total angular momentum I. When the effective nucleon-nucleus interaction V can be 
expanded in irreducible tensors one finds that (ref. 2) 



= (21 + 1 ) 


2Jj + 1 


2J 2 + 1 


1/2 


LSJ 


X 


U(j 2 IJJ-p: ^1^ ( a 2 J 2 ll ^LSj( r ’ (®) 


where (a^J^) and (ag^) refer to the initial and final nuclear states, respectively, and 
j^ and j 2 to the initial and final projectile states. The reduced matrix elements of 

V L gj(r, A) can be obtained, when the nuclear wave functions are of the projected HF 
type, by the methods outlined in reference 8. Matrix elements of V diagonal in the 
channel index c are parameterized with the standard optical model potential (ref. 8) 
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U opt< r > = V c< r > - V 0 *>V< r) + 4iW S a S “ »S^ + V LS 



g • l d 
r dr 


P LS^ 


(10) 


The nondiagonal elements are calculated with a nucleon-nucleus interaction that is a 
sum of potentials representing the nucleon-nucleon interaction. These two-body poten- 
tials have the simple form (ref. 1) 


v(r) = -52 exp 



6 p se) 


(ID 


where P-pg and Pgg are projection operators for the triplet-even and singlet-even 
states, respectively. We recognize that even for the nondiagonal matrix elements, a 
renormalized interaction should be used; however, we feel that the preliminary nature 
of this investigation does not warrant the additional complexity. In addition we have 
neglected the knock-out exchange amplitudes (refs. 10 and 11) since we feel, in many 
instances, their importance has been overestimated. 

Nuclear structure information enters the scattering calculation through structure 
amplitudes characterizing an elementary transition between nuclear states. These 
quantities are simply related to the Sj's which were defined in equation (6) and which 
express the likelihood that a transition will proceed via a scattering through the config- 
uration (a, b). The structure factors have been calculated for HF model spaces com- 
posed of six, 10, and 15 orbitals in order to study the way in which truncation of the 
model space affects inelastic scattering. These spaces are defined in table I. 


TABLE I. - DESCRIPTION OF BASIS SPACES TO 
BE USED IN HARTREE-FOCK CALCULATIONS 


Space 


Basis 

( ls l/ 2 , 2Sjy 2 , ld 3/ / 2 , ld 5 / 2 ) dPi/2- 1 P3/2 ) 

( ls l/2’ 2 s 1/2’ ld 3/2’ ld 5/2 ) 

< lp l/2’ lp 3/2’ 2p l/2’ 2p 3/2’ lf 5/2- lf 7/2 ) 

(Is 1/2 ’ 2Sl 


T/2- ld 3/2' ld 5/2’ 3 s 1/2' 2d 3/2’ 
2d 5/2’ lg 7/2' lg 9/2^ 

( lp l/2 ’ lp 3/2 ’ 2p l/2 ’ 2p 3/2' lf 5/2 ’ lf 7/2 ) 
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RESULTS WITH SIMPLE HARTREE-FOCK MODEL 

The simple Hartree-Fock model has received a great deal of attention in the last 
few years. The prediction of many gross nuclear properties of light nuclei, based on 
the use of this model has been quite successful (refs. 1, 7, and 12). However, use of 
even the largest available model spaces in such studies seems to fall short in predict- 
ing the magnitudes of various transition rates, and intrinsic quadrupole moments. 
Bassichis et al (ref. 12) have found that model spaces containing as many as 21 basis 
states do not yield intrinsic quadrupole moments which have converged, and Ford, 
Braley, and Bar-Touv (ref. 13) have found similar results for projected E2 rates 
using model spaces containing up to 15 orbitals (refs. 12 and 13). In spite of these tem- 
porary short comings, the value of the simple HF approximation as a research tool is 

20 24 

quite apparent and has been used to describe the nuclei of interest here ( Ne and Mg). 

20 24 

Wave functions for Ne and Mg, which are used to generate nuclear form 
factors, are obtained from reference 10. The relevant model spaces are defined in 
table I. A brief review of some of the structure results from reference 10 appear in 
tables 11(a) and (b). (The projected energy spectra may be found in figs. 1 and 2. ) The 
oscillator parameters (b Vhu>/ m) used in these cases were chosen to reproduce the 
experimental RMS radii since the two-body force used in the structure calculations was 
nonsaturating and the size of the basis was not sufficiently large to provide the correct 
long range behavior of the wave function (ref. 13). In both cases the use of a more 
realistic two-body force would probably result in smaller oscillator lengths to repro- 


TABLE II - COMPARISON OF PREDICTED AND EXPERIMENTAL 
RADII AND E2 RATES 
(a) Neon-20 



<r 2 > l / 2 . 

B(E2; 0 + - 2 + ), 

B(E2 ; 2 + - 4 + ), 

B(E2; 4 + - 6 + ), 


fm 

e 2 • fm 4 

2 , 4 

e ■ fm 

e 2 • fm 4 

Predicted 

2. 78 

204 

103 

90 

Experiment 

2.79 

286±15 

89±9 

129±13 


(b) Magnesium-24 



(R 2 ) 1 ' /2 , 

B(E2; 0 + - 2 + ), 

B(E2; 2 + - 4 + ), 


fm 

e 2 ■ fm 4 

e 2 ■ fm 4 

Predicted 

3.04 

374 

191 

Experiment 

3. 02±0. 03 

436±46 

133±13 
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PHF EXP'T 


Figure 1. - Experimental and projected 
Hartree-Fock energy spectrum for 
neon-20 with b = 1.93 fm. 



PHF EXP'T 


Figure 2. - Experimental and projected 
Hartree-Fock energy spectrum for 
magnesium-24 with b=2.09fm. 


duce the same radii. It is difficult to predict how this might affect the transitions rates 
and other properties that are known to be sensitive to the long range behavior of the 
wave function. 

Coupled channel and DWBA calculations were carried out for 24. 5-MeV proton 
scattering from ^Ne. The o + , 2 + , 4 + , and 6 + members of the ground-state rotational 
band were included in the coupled equations. The optical parameters used in the CC 
calculation were adjusted so the CC elastic scattering gave the best fit to the experi- 
mental elastic data, whereas the DWBA optical parameters were adjusted to yield the 
same elastic scattering as the CC result (see fig. 3). The CC calculation is thus being 
used as a standard against which the DWBA calculation may be compared. The optical 
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Figure 3. - Comparison of elastic scattering of 24. 5-MeV protons from neon-20 predicted by coupled- 
channel calculation and standard optical model calculation with DWBA parameters described in the 
text. 


TABLE in - COMPARISON OF CC AND DWBA OPTICAL 
MODEL PARAMETERS 



V o’ 

W s’ 

V LS’ 

V 

r s’ 

r LS’ 

V 

a s’ 

a LS’ 


MeV 

MeV 

MeV 

fm 

fm 

fm 

fm 

fm 

fm 

Neon-20: 










CC 

55. 4 

7. 05 

3. 58 

1. 05 

1.265 

0. 96 

0. 73 

0. 61 

0. 33 

DWBA 

54. 5 

8. 12 

3. 58 

1. 073 

1.251 

.96 

.713 

. 61 

.33 

Magnesium -24: 










CC 

49.1 

6. 56 

5. 29 

1. 174 

1. 19 

1. 06 

0. 736 

0. 562 

0.546 

DWBA 

46. 1 

8.30 

5. 29 

1. 241 

1. 163 

1. 06 

.665 

. 552 

. 546 


parameters are listed in table HI; the main difference, as one would expect, is in the 
imaginary surface derivative potential (W g ), which has changed by 13 percent in going 
to the DWBA limit. 

A comparison between CC and DWBA calculations in the largest model space 
(fig. 4) shows that the cross section for excitation of the 2 + level is hardly changed in 
the DWBA limit. This is to be expected since the 2 + level is predominantly coupled to 
the ground state. If one compares cross sections the 4 + levels, however, the CC result 
is enhanced by a factor of about 1. 7, and the angular distribution is somewhat different, 
particularly at larger angles. These changes are due to the double excitation through 
the 2 + levels that cannot be accounted for in a DWBA calculation. Since the 6 + level 
is fed by both the 2 + and the 4 + levels, its cross section shows an even larger deviation 
between DWBA and CC predictions. The discrepancy in magnitude between theory and 
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Figure 5. - Coupled channel results compared with experiment for 
scattering of 24. 5-MeV protons from neon-20. 


experiment (fig. 5) is not surprising since a bare nucleon-nucleon force is used in the 
scattering problem. 

We have also studied how the cross sections are affected by truncation of the 
nuclear model space. Coupled channel and DWBA calculations were carried out with 
model spaces containing six, 10, and 15 basis functions, as described in table I (see also 
fig. 6). For the sake of simplicity these calculations were made with the same set of op- 
tical parameters which appear in table III. As one can see from figure 7, the major en- 
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Cross section, mb/sr 



Figure 6. - Coupled-channel prediction of elastic scattering of 24. 5-MeV protons from 
neon-20 compared with experiment. The nuclear model spaces used in the calculation 
contain six states and 15 states as indicated. 


Cross section, mb/sr 



Angle, 6 cm , deg 


Figure 7. - Comparison of coupled-channel and DWBA predictions with experi- 
ment for excitation of the 2 + (Q * -1. 63 MeV) state in neon-20 by 24. 5-MeV 
protons. The effect of varying the size of the model space is also indicated. 



hancement of the 2 + level cross section comes from the inclusion of the 2p - If shell. No 
such saturation seems to be in sight for the 4 + level (fig. 8). In fact, going from the 10- 
state to the 15-state basis increased the 4 + cross section by almost a factor of two. 

(The plot of the 10-state result is not included because it would become difficult to read 
the graph). This is probably due to the inclusion of the lg- 7/2 an( * ^9/2 s ^ es that 
carry four units of orbital angular momentum. The magnitude of the cross section for 
the 6 + level (fig. 9) is insignificant until one includes at least fifteen basis functions in 
the model space. 

20 

Similar calculations as described for Ne were carried out for the scattering of 
17.5-MeV protons from ^Mg. The 0 + , 2 + , and 4 + members of the ground-state rota- 
tional band were included in the coupled equations. As before, the optical parameters 
used in the CC calculation were adjusted to yield a best fit to the available experimental 
elastic data (ref. 6), as shown in figure 10. The DWBA parameters were again chosen 



Figure 10. - Comparison of elastic scattering predicted by coupled-channel calculation and standard 
optical model calculation with DWBA parameters described in the text. 


to yield the same elastic scattering as the CC result; see table IE. 

The DWBA and CC results for inelastic scattering are compared in figure 11. As 
20 + 

in the case of Ne, the 2 cross sections are seen to be quite similar; however, the 
CC 4 + cross section is now enhanced by almost a factor of 10 over the DWBA result. 
Thus, it is likely that a considerable amount of the flux to the 4 level in Mg is due to 
a double excitation process. The CC results are compared with experiment in figure 12. 
As discussed earlier, the discrepancy in magnitude between theory and experiment is 
not surprising. This will be discussed in greater detail in the following section. 
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Cross section, mb/sr 



Angie, 0 cm , deg 


Figure 11. - Comparison of DWBA and coupled-channel predicted cross 
sections for the scattering of 17.5-MeV protons from magnesium-24. 



Angle, e cm , deg 

Figure 12. - Coupled-channel results compared with experiment for 
scattering of 17. 5-MeV protons from magnesium-24. 




INSIGHTS FROM THE MACROSCOPIC MODEL 


The CC-predicted cross sections for excitation of the 2 + levels in both nuclei have 
good angular shapes when compared with experiment, and as one would expect the CC 
and DWBA cross sections for these levels are in close agreement. However, such is 
not the case for the 4 + levels in these nuclei, as was noted earlier. We have found that 
considerable insight into this difference may be gained by examining the results of CC 
calculations made with the macroscopic model (ref. 6). These predictions, obtained 

using the collective model, appear in figures 13 and 14. 

20 

Examination of the Ne results reveals that when only quadrupole shape correla- 
tions are included, that is, /3^ = 0, the predicted cross sections do not compare favor- 
ably with experiment. Moreover, the quality of the shape and magnitude deteriorates 
as the excitation energy increases. 

The next highest order in the deformation involves the hexadecapole shape param- 
eter Some cross sections calculated using both quadrupole and hexadecapole terms 

^Netp, p'l^Ne* 

E p = 24. 5 MeV 


Experiment 

po = 0. 47 p 4 = 0 



Angle, 0 cm , deg 


Figure 13. - Coupled-channel calculation for Neon-20 
based on the macroscopic model compared with ex- 
periment. (See ref. 11.) 
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Figure 14. - Coupled-channel calculation for 
magnesium-24 based on macroscopic model 
compared with experiment. (See ref. 11.) 


also appear in figure 13. The excitation of the 2 + level is not affected significantly, but 

it is evident that the 4 + cross section cannot be fit unless the deformation parameter /3 4 

2 ’ 
is given a positive value of the order of This implies that the indirect excitation 

(0 + - 2 + - 4 + ) is roughly as effective as the direct excitation (0 + - 4 + ) in contributing to 

the 4 + cross section. The DWBA calculation includes only the latter, and so may be 

expected to underestimate the CC value by about a factor of 2. This supposition, a 

direct inference from macroscopic model results, is borne out in our microscopic 

model calculations by the 4 + cross sections displayed in figure 4. 

Another piece of evidence may be obtained by using the fact that, judging from the 
macroscopic calculations, excitations of the 6 + level are the most sensitive to hexadeca- 
pole deformation. A comparison of the collective model 6 + cross sections shown in fig- 
ure 13 to the Hartree-Fock model result shown in figure 4 suggests that the microscopic 
wave function does not have enough hexadecapole deformation and instead corresponds 
roughly to ( 3 ^ = 0 . This in turn suggests that improvement in the microscopic wave 
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function could be obtained if a means were found to enhance the hexadecapole deforma- 
tion. 

24 

A somewhat different situation prevails in the case of Mg. Here, the best macro- 
scopic model fits to the data (fig. 14) are achieved for deformation parameter values 
such that [/3^|« 02 . Consequently one would expect the DWBA calculation of the 4 + 
cross section, which treats only the direct excitation, to considerably underestimate the 
CC result. This indeed occurs, as may be seen in figure 11. Once again the HF wave 
function has demonstrated characteristics known to be necessary from the macroscopic 
studies. 

We have seen that calculations of scattering cross sections performed using pro- 
jected HF nuclear wave functions exhibit the same general behavior as calculations using 
the macroscopic model. This indicates that the moments, or shape, of the nucleus 
must be quite well described by the HF model. One last confirmation of this is provided 

by table IV, in which the nuclear moments in the intrinsic state are listed; the small- 
24 

ness of the Mg hexadecapole moment is striking. 


TABLE IV. - MASS QUADRUPOLE 
(Q 2 ) AND HEXADECAPOLE 
(Q 4 ) MOMENTS 



Q 2 , fm 2 

Q 4 , fm 4 

20 Ne 

83. 5 

304. 9 

24 Mg 

114. 5 

128. 6 


CONCLUSIONS 

The primary object of this investigation has been to make a detailed comparison of 
DWBA and CC predicted cross sections and to examine, explicitly, how scattering cal- 
culations are affected by truncating the bound-state model space. The reliability of 
such a study ultimately depends on comparison with experiment. It is at this point that 
the decision must be made whether one is testing the nuclear model, the scattering 
mechanism, or the two-body force. In contrast to particle -transfer reactions, the 
mechanism for inelastic scattering of nucleons is thought to be well understood; further- 
more, the Glendenning-Veneroni two-body interaction has been tested in earlier studies 
(ref. 14) and has yielded accurate angular distributions (although it should be mentioned 
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that Glendenning and Veneroni found it necessary to multiply their cross sections by 
factors from three to eight in order to fit the data). Thus, it was our original intention 
to use inelastic proton scattering as a means of testing the reliability of the best avail- 
able Ne and Mg nuclear wave functions and to determine at the same time those 

situations in which the DWBA may be used instead of the considerably more involved CC. 

. 20 

The DWBA appears to be adequate for excitation of the 2 level in both Ne and 

24 Mg, but cross sections for the higher lying states are altered significantly by the CC 

calculation. In particular, the DWBA and CC predictions for the 4 + cross section in 
24 

Mg are very different, which can be understood in terms of the deformation param- 
eters used in the macroscopic model. 

j_ 24 

A very similar difference in cross section for the 4 level in Mg was also obtained 
using the projected Hartree-Fock nuclear wave function, which provides gratifying evi- 
dence of the reliability of this microscopic nuclear model. After all, analysis of nuclear 
reactions by means of a macroscopic model essentially yields information about gross 
features of the nucleus, namely, its quadrupole and hexadecapole moments, as measured 
by the deformation parameters /3g and (}.. It now seems clear that these same general 
features can also be reproduced by the much more complicated many-nucleon wave 
function, which is determined solely from general symmetry constraints and the require- 
ment of energy minimization. 

There remain some areas of uncertainty, however. Including the active core (which 
was left out of Glendenning and Veneroni 's calculation), enlarging the model space, and 
using the coupled channel method have all helped to bring the cross sections closer to 
agreement with experiment - but not enough. It seems unlikely that further enlargement 
of the bound-state model space will solve the problem. One is therefore led, somewhat 
reluctantly, to question the validity of the two-body force being used, and perhaps even 
the inelastic scattering mechanism. It is well known that the two-body force is stronger 
within nuclear matter than it is in vacuum, and it would not be unreasonable to expect 
this renormalized force to be responsible for the nuclear transition. It is also possible 
that effects of antisymmetrization between the projectile and the target nucleons, usually 
considered small and therefore ignored as part of the excitation mechanism, may in 
fact play a substantial role. Both these possibilities are currently under investigation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 15, 1971, 

129-02. 
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